    
dir1 = ('./figure/');
colormap jet
    set(gcf,'position',[500,500,1000,200])
    pcolor(yy(2:end-1),xx(1:end),squeeze(T0i(1:end,2:end-1)+Ti(1:end,2:end-1)-1)); colorbar('Fontsize',15); shading interp;
    title(['T,t=' 24000],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
%     caxis([-2 2])
    drawnow
    print(gcf,'-dpng',[dir1,sprintf('TT%d',56)]);
    close